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ABSTRACT 

Scattering and absorption properties at optical and ultraviolet wavelengths are cal- 
culated for an interstellar dust model consisting of carbonaceous grains and amorphous 
silicate grains. Polarization as a function of scattering angle is calculated for selected 
wavelengths from the infrared to the vacuum ultraviolet. 

The widely-used Henyey-Greenstein phase function provides a good approximation 
for the scattering phase function at wavelengths between ~0.4 and 1/im, but fails to fit 
the calculated phase functions at shorter and longer wavelengths. A new analytic phase 
function is presented. It is exact at long wavelengths, and provides a good fit to the 
numerically-calculated phase function for A > 0.27/iin. 

Observational determinations of the scattering albedo and (cos 9) show considerable 
disagreement, especially in the ultraviolet. Possible reasons for this are discussed. 

Subject headings: dust, extinction - polarization - scattering - ultraviolet: ISM 

1. Introduction 

Interstellar grains scatter electromagnetic radiation. Reflection nebulosities are conspicuous at 
optical and UV wavelengths when dust is brightly illuminated by a nearby star. Dust clouds which 
are not unusually close to a star are illuminated by the general interstellar radiation field. Finally, 
the starlight scattered by dust in the diffuse interstellar medium constitutes the so-called "diffuse 
galactic light". Observations of reflection nebulae, dust clouds, and the diffuse galactic light, provide 
a means of determining the scattering properties of interstellar grains, thereby testing models for 
interstellar dust. 

Photoionization and photodissociation of molecules play a major role in interstellar chemistry, 
and the chemical structure of molecular clouds is therefore directly linked to the ability of ultraviolet 
starlight to penetrate into these dusty regions. Knowledge of dust scattering properties in the 
ultraviolet is therefore required for realistic modeling of interstellar clouds. 

The nature of interstellar grains remains uncertain (see Draine 2003a, and references therein). 
This paper will examine the scattering properties for a grain model consisting of two separate grain 
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populations - carbonaceous grains and silicate grains. With the grains approximated by homoge- 
neous spheres with the size distributions found by Weingartner & Draine (2001; hereafter WD01), 
this grain model is consistent with the observed interstellar extinction, the observed infrared emis- 
sion from interstellar dust (Li & Draine 2001, 2002), and the X-ray scattering halo observed around 
Nova Cygni 1992 (Draine & Tan 2003). The carbonaceous grains are assumed to be primarily car- 
bon when the grains are large, but to extend down to very small sizes with the smallest grains being 
individual polycyclic aromatic hydrocarbon molecules. The scattering is dominated by grains with 
radii a > 100 A, containing > 10 6 atoms; carbonaceous grains in this size range are modeled using 
the optical properties of graphite. 

The primary objective of this paper is to calculate the scattering and extinction properties 
of this dust model at infrared, optical, and ultraviolet wavelengths, and to make these results 
available for use in radiative transfer calculations and for comparison with observations. The X-ray 
scattering and absorption properties of this grain model are the subject of Paper II (Draine 2003b). 

The adopted dielectric functions are presented in §2. The scattering properties of interstellar 
dust at optical and ultraviolet energies, as calculated for the carbonaceous-silicate grain model, are 
presented in §3. We show the scattering phase function at selected wavelengths from the SDSS z 
band (A = 8930 A) to the vacuum ultraviolet (A = 1820 A). Scattering properties are calculated 
for Milky Way dust with Ry = 3.1, and also for models for dust in the LMC and SMC. In §3.2 we 
show the degree of polarization as a function of scattering angle for selected wavelengths. 

The Henyey-Greenstein phase function has often been used to approximate the anisotropic 
scattering properties of interstellar dust. In §4.1 we show that the Henyey-Greenstein function has 
r.m.s. error < 10% for 0.47/im < A < 0.94/um, but has larger errors outside this range. We present 
a new analytic phase function (equation 5) with a wider range of applicability, with r.m.s. error 
< 10% for A > 0.27/xm. 

In §5 we collect observational determinations of the scattering albedo and (cos 6) for dust in 
reflection nebulae, in dense clouds, and in the diffuse interstellar medium. Discrepancies among 
these determinations are noted, and possible reasons for this are discussed. 

The principal results are summarized in §6 

2. Dielectric Function 

As discussed by WD01, the grain population is assumed to include very small grains with the 
optical properties of polycyclic aromatic hydrocarbon molecules (PAHs), plus larger grains which 
are approximated as carbonaceous or silicate spheres. The PAHs produce negligible scattering. 

From the observed 3.4/um C-H stretch feature, Pendleton & Allamandola (2002) estimate that 
~85% of the C is aromatic, and ~15% is aliphatic (chainlike). The graphite dielectric function 
will be used to approximate the optical and ultraviolet response of interstellar carbonaceous grain 
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material. Scattering and absorption by the carbonaceous spheres is calculated using the dielectric 
tensor of graphite, using the usual "1/3-2/3 approximation" (Draine & Malhotra 1993). 

The dielectric functions used here are taken from Paper II, which constructs self-consistent 
dielectric functions extending from microwave to X-ray energies, including realistic structure near 
X-ray absorption edges. The adopted dielectric functions for graphite and "astronomical silicate" 
are shown in Figures 1 and 2. These dielectric functions are close to the dielectric functions obtained 
previously by Draine & Lee (1984), although differing in detail. 
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Fig. 1. — Dielectric function for graphite (Draine 2003b). 

3. IR-Optical-UV Scattering by Interstellar Dust 

3.1. Angular Distribution of Scattered Light 

Weingartner & Draine (2001; hereafter WD01) obtained size distributions of spherical carbona- 
ceous and silicate grains which reproduce the observed extinction curve both in the local Milky 
Way and in the Large and Small Magellanic Clouds. Here we calculate the scattering properties of 
these dust mixtures. 
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Fig. 2. — Dielectric function adopted for amorphous MgFeSiCU (Draine 2003b). 

The scattering properties of a particular dust mixture X are characterized by the differential 
scattering cross section per H nucleon 




where (drij /da)da is the number of grains of type j per H nucleon with radii in (a, a + da), and 
(dC sca /dCl)j ta; \ is the differential scattering cross section for grain type j, radius a, at wavelength 
A, for scattering angle 6, for a grain illuminated by unpolarized light. The grains are assumed to 
be spherical, and the differential scattering cross sections are calculated using Mie theory (Bohren 
& Huffman 1984), using the code developed by Wiscombe (1980, 1996). 

In Figure 3 we show the differential scattering cross section per H nucleon for Ry = 3.1 Milky 
Way dust, at the central wavelengths of SDSS z (8930A), i (7480A), r (6165A), g (4685A), and 
u (3550A), Cousins I (8020A) and R (6492A), V (5470A), and the F250W (2696A), F220W (2220A), 
and F25CN182 (1820A) filters for the Space Telescope Imaging Spectrograph (STIS). The scattering 
becomes stronger and more forward-throwing at shorter wavelengths. 

To see the sensitivity to variations in the dust mixture, Fig. 3 also shows da SCSL /d£l calculated 
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Fig. 3. — d<7 aca /dr2 at selected wavelengths A, as a function of the scattering angle 6, for WD01 models for Milky 
Way dust with Ry = 3.1, LMC average dust, and SMC bar dust. Curves are labelled by wavelength A. 



for the WD01 dust mixtures for the "average LMC" and the SMC bar. The calculated scattering 
closely resembles the Milky Way scattering, but with an overall reduction of about a factor of 4 for 
the LMC, and a factor of 8 for the SMC, in line with the overall dust and metal abundance relative 
to the Milky Way. 1 

Figure 4 shows the wavelength dependence of <7 ext (the total extinction cross section per H), 
the scattering albedo <7 sca /<7 cx t, and the first and second moments 2 (cos 6) and (cos 2 8) for the 



^e/H is ~ 30% of solar in the LMC and ~ 14% of solar in the SMC (Dufour 1984; Kurt & Dufour 1998). In the 
LMC, E(B - V)/Nh is ~ 24% of the local Milky Way value (Koorneef 1982; Fitzpatrick 1986) while in the SMC it 
is only ~ 13% of the local value (Martin et al. 1989). 

2 The second moment (cos 2 ) will be used in the new analytic phase function proposed in §4.2. 
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Fig. 4. — Extinction cross section per H, albedo, (cos), and (cos 2 ) calculated for WD01 models for Milky Way dust 
with Rv = 3.1, 4.0, and 5.5, average LMC dust, and dust in the SMC bar. 

scattered radiation. Six different grain models are shown, fitted to different observed extinction 
curves (see WD01 for details). There are considerable differences in albedo, (cos), and (cos 2 ) among 
the models. In particular, the SMC bar model shows a high albedo near 4.6//m _1 , in contrast to 
the other five models which have a local minimum in the albedo at this frequency. This is because 
the SMC bar model differs from the other 5 models in lacking PAHs and a < 0.02/im graphite 
grains, as these are not allowed by the observed absence of a 4.6/im -1 extinction "bump". 
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3.2. Polarization of Scattered Light 

Even when a grain is illuminated by unpolarized light, the scattered radiation is generally 
polarized. The degree of polarization depends upon both the scattering angle and the wavelength 
of the radiation. The fractional polarization P = (I± — (I± +I\\), where I± , Jy are the intensities 
of scattered light in polarization modes perpendicular or parallel to the scattering plane. 

The polarization P of the scattered light is shown in Fig. 5 as a function of scattering angle 9 
for Milky Way dust with Ry = 3.1, at 11 different wavelengths. Rayleigh scattering would have P = 
(1 — cos 2 #)/(l+cos 2 9), with P = 1 for 9 = 90°. At long wavelengths, the polarization has a distinct 
peak near ~ 90°, but even for A ~ 0.9/im the peak polarization is only ~ 0.45. As the wavelength 
is reduced, the peak polarization declines. For A < 0.6/um, the polarization becomes negative at 
large scattering angles (see Figure 5), with large negative polarizations in the 120-150° region for 
0.2//m < A < 0.4/zm. Note, however, that this negative polarization occurs for scattering angles 
where the scattering is very weak (see Figure 3) and therefore could be masked by scattering at 
other points on the sightline where the scattering contributes a positive polarization. Observations 
of the predicted negative polarization will probably require simple scattering geometries, such as 
dust in a thin disk, illuminated by a single source. 

Similar results are found for the LMC and SMC dust mixtures - see Fig. 5. Note the very large 
negative polarizations found for the LMC mixture for 0.2/xm < A < 0.55/xm. From the variation of 
the ultraviolet polarization signature between the different grain size distributions in Figure 5 it is 
apparent that the ultraviolet polarization is sensitive to the details of the grain size distribution. 



4. Analytic Approximations for the Phase Function 

4.1. Henyey-Greenstein Phase Function 

Scattering of unpolarized incident light by a dust mixture at wavelength A is characterized by 
the total scattering cross section per H nucleon, cr sca (A), and a "phase function" 

characterizing the angular distribution of the scattered light, where da sc& /dVL is the differential 
scattering cross section, at scattering angle 9, for unpolarized incident light. The definition (2) 
implies the normalization J dQ <3? = 1. Isotropic scattering would have $ = 1/47T. The first 
moment 

(cos 9) = J dQcos9 $(0,A) (3) 
of the phase function is a measure of the asymmetry between forward and backward scattering. 
Henyey & Greenstein (1941) proposed an analytic function to model anisotropic scattering for 
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Fig. 5. — Degree of polarization as a function of scattering angle 8, for scattering by Milky Way dust with Rv = 3.1, 
LMC average dust, and SMC bar dust. Curves are labelled by wavelength A. 



dust grain mixtures: 



4vr (1 + g 2 - 2g cosO) 3 / 2 



9 = (cos 6) . (4) 



With a single parameter g, this is a convenient analytic form that has been widely used to represent 
dust scattering properties in radiative models of dusty regions. 

By construction, 4>o{6) has the correct first moment (cos 6) = f dQcosO (j>o(0) = g, but of 
course 4>o does not perfectly reproduce the angular dependence of the real phase function As 
will be seen below, <j)Q is a poor approximation at both A < 0.4/im and A > 1/xm. 
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4.2. 



A New Phase Function 



At long wavelengths, Rayleigh scattering prevails, with (cos9) — > and $ — > (3/167r)(l + 
cos 2 #). When (cos 9} — > 0, the Henyey-Greenstein phase function </>o — > 1/4-7T; </>o is 33% low at 
6* = 0, 7T and 33% high at 9 = vr/2. 

Consider the phase function 



with two adjustable parameters, a and g. For a = 0, equation (5) reduces to the Henyey-Greenstein 
phase function <j)Q (equation 4) . For g = and a = 1 this reduces to the phase function for Rayleigh 
scattering. For a = 1 this corresponds to the phase function proposed by Cornette 8z Shanks (1992). 
Analytic results for this phase function are given in Appendices A - C. 

To determine the parameters a and g we can require that <p a (9) have the same first and second 
moments (cos9) and (cos 2 9) as <3>. The parameters g and a are then given by equations (B4-B6). 
However, although the resulting phase function has correct first and second moments, the fit is 
poor when the dust is strongly forward-scattering ((cos#) > 0.6). 

We find that an improved fit is obtained if we obtain g and a from equations (B4-B6) only 
when the resulting a < 1; for values of (cos 9} and (cos 2 9} for which equations (B4-B6) lead to 
a > 1, we instead set a = 1 and obtain g from equation (C2). We will refer to this new analytic 
phase function as 4> a <\. 

Figure 6 shows the scattering phase function 3>(0) for the WD01 Milky Way dust model, 
together with the Henyey-Greenstein phase function and our new phase function </> Q <i, at wave- 
lengths A = 1.2/im, 0.80/um, and 0.40/im. For these three wavelengths, the new phase function 
provides an improved fit to the actual phase function $. For A = 0.40/um, we also show <3? ap- 
proximated by a sum over Legendre polynomials Pi up to I = 8. This 9 term expansion does not 
reproduce $ as well as the new phase function 4> a <i. 

At ultraviolet wavelengths the grains become more forward-throwing, and the analytic phase 
functions <j>o and 4> a <i no longer provide a good fit, as seen in Figure 7 for A = 0.20^m and O.lO/um. 
At A = 0.20/mi, for example, both the Henyey-Greenstein phase function </>o and the new phase 
function (j) a <i underestimate the forward scattering intensity by a factor of ~2 for cos 9 > 0.98 



To quantify the error associated with using an analytic phase function (f)(9) to approximate an 
actual phase function &(9), we define the r.m.s. relative error 



M0) = -r- 



1 1 - g 2 1 + a cos 2 9 



(5) 



4vr [l + a(l + 2 5 2 )/3j (l + g 2 - 2g cos 9 fl 2 



{9 < 10°). 




(6) 
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Fig. 6. — Scattering phase function for the WD01 Milky Way dust model at three wavelengths (solid lines) compared 
to the Henyey-Greenstein phase function (f>o and our new phase function <f> a <i. For A = OAfim a 9 term Legendre 
polynomial representation is also shown. 



and r.m.s. absolute error 



^■abs — 



4-7T 



1/2 



(7) 



where, of course, (<&) = 1/4tt. The r.m.s. relative error h Te \ would appear to be the best measure 
of the overall quality of fit. h Te \ and /i a b s will differ substantially only when the phase function is 
very asymmetric; large fractional errors in directions where the scattering is weak then make only a 
small contribution to /i a bs> while modest fractional errors in directions where the scattering is very 
strong will make large contributions to /i a bs- 
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Fig. 7. — Same as Figure 6, but for A = 0.2/im and 0.1/^m. Legendre polynomial representations of <fr(8) are shown 
with 16 and 51 terms. 



In the limit of Rayleigh scattering (A 3> a), &(0) — ► | (l + cos 2 6>), and 4>o has errors 



^abs 
^rel 



1/V20 » 0.22 
(13 - 4vr) 1 / 2 /3 



0.22 



(8) 
(9) 



Figure 8 compares the error of different analytic approximations to the phase function, using 
as an example <&(#, A) for the WD01 model for Ry = 3.1 Milky Way dust. The errors h Te \ and /i a b s 
are shown in Figure 8 for the Henyey-Greenstein function 0q, the Cornette-Shanks phase function 
(pi, and our new phase function <j) a <\. At long wavelengths A > 3fim, the Henyey-Greenstein 
approximation fa has h ie \ « h a ^ s 0.22, as expected from eq. (8,9). Between ~ 0.5 and ~ lfim, 
the relative errors for <pQ are modest, h TC \ < 10%, so (f>o provides a good approximation to the actual 
scattering properties. However, cpQ has h TC \ > 10% for A < 0.5//m, rising to shorter wavelengths. 
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Fig. 8. — Normalized relative error h le \ (eq. 6) and absolute error /i a b s (eq. 7) for the analytic phase functions <j>o, 
4>i, and 4> a <i applied to the WD01 grain model for MW dust with Ry = 3.1. 

Figure 8 shows the fractional errors /i re i and hghs f° r the new phase function <j> a <\. As expected, 
the new phase function provides an excellent fit at long wavelengths, with fractional errors < 3% 
for A > l//m. As the wavelength is decreased, the new phase function (j) a <i (as measured by h re \) 
remains preferable to the Henyey-Greenstein phase function for A > 0.19/zm. <p a <i has h Ie \ < 10% 
for A > 0.27/xm. 

In the ultraviolet the phase function becomes very strongly forward-scattering, and all of the 
the analytic approximations have significant errors. More complicated parameterizations - such 
as the use of multiple Henyey-Greenstein phase functions, as by Witt (1977) and Hong (1985) - 
may be considered. Alternatively, a given phase function <£(#) can be approximated by summing 
over Legendre polynomials, although the sum needs to include a large number of terms to provide 
a good approximation - 16 terms are required for A = 0.20/xm, and even 51 terms are insufficient 
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for A = 0.10/um. At short wavelengths it may be best to simply tabulate 3>(0) and interpolate as 
required. 

5. Discussion 

The scattering properties of interstellar dust have been determined observationally by com- 
paring the observed surface brightness of reflection nebulae with model nebulae computed with 
different dust scattering properties, selecting the model which provides the best match to the ob- 
servations. The usual approach has been to try to determine only two grain properties - the albedo 
and (cos 6) - by finding radiative transfer models which appear to be consistent with the observed 
surface brightness of scattered light. 

It has been customary to assume that the scattering phase function <J>(0) can be approximated 
by the Henyey-Greenstein function 4>o{9), with g = (cos 9) 

This approach has been used with (1) observations of high surface brightness reflection nebulae 
(e.g., NGC 2023, NGC 7023, IC 435) illuminated primarily by a single star, (2) observations of 
individual dust clouds illuminated by ambient starlight, and (3) observations of the much fainter 
diffuse galactic light (the entire Galaxy as a reflection nebulosity). Figures 9 and 10 show the results 
so obtained from a number of independent studies, at wavelengths from the optical to vacuum 
ultraviolet. Also shown are the albedo and asymmetry factor calculated for WD01 grain models 
representing average Milky Way dust with Ry = Ay/E(B — V) = 3.1, and dust from denser 
regions with Ry = 4.0 and 5.5. While there are differences among the different observational 
determinations, particularly in the ultraviolet, many of the observational results are in general 
agreement with the albedo and asymmetry factor calculated for the WD01 grain model. 

For A -1 < 4 / [/m~ 1 , the observational studies appear to be in general agreement with one 
another, but at shorter wavelengths the observational results are sometimes in conflict. For example, 
from observations of the diffuse galactic light at ~ 6.25//m _1 (A = 0.16/xm) Hurwitz et al. (1991) 
found (a,g) = (0.185±.055,0.2±0.2), while Witt et al. (1997) found (a,g) = (0.45±0.05, 0.68±0.10). 

Aside from genuine regional variations in dust properties, there are many factors which could 
contribute to such discrepancies. Foremost may simply be the difficulty of observations in the vac- 
uum ultraviolet, followed by uncertainties concerning the nebular geometry and the illuminating 
radiation. However, given that use of the Henyey-Greenstein phase function <pQ introduces an abso- 
lute error /i a b s > 70% in the scattering phase function at A < 0.16^m (see Fig. 8), it seems possible 
that reliance on the Henyey-Greenstein phase function may contribute to the disagreements among 
different observational determinations of the albedo and scattering asymmetry in the ultraviolet. 

Unfortunately, the new phase function (j) a <i is inaccurate for A < 0.27/um. For A -1 > 4/im" 1 , 
radiative transfer models should use a more accurate representation of the phase function <&. As 
discussed above, one possibility is expansion in Legendre polynomials, although many terms are 
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Fig. 9. — Scattering albedo as a function of frequency for different MW dust mixtures, together with observational 
estimates of the dust albedo in various regions. 



required. Alternatively, it may be simplest to use tabulated values of <& with interpolation. "Inver- 
sion" of the observational data to infer the grain properties may not be feasible; the best approach 
may be to assume the dust properties [i.e., (Text (A), albedo(A), and $(6*, A)] based on a grain model, 
and then seek a dust spatial distribution which maximizes consistency with the observed surface 
brightness and whatever other constraints are available. Failure to find agreement would suggest 
that the grain model may be incorrect. 

Applying this approach to reflection nebulae is complicated by the possibility of clumpy dust 
distributions, which Mathis et al. (2002) show can produce a range of ratios of nebular surface 
brightness to unscattered stellar flux, depending on the detailed dust clump distribution and the 
viewing direction. Diffuse clouds (i.e., the diffuse galactic light) and externally-illuminated clouds 
may prove easier to interpret than bright reflection nebulae, which generally have a bright star 
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Fig. 10. — Scattering asymmetry factor g — (cos 9) as a function of frequency for different MW dust mixtures. Also 
shown are observational estimates for g in various regions. 



embedded in the dust distribution. 



6. Summary 

The following are the principal results of this work: 

1. Differential scattering cross sections have been calculated at selected wavelengths from in- 
frared to ultraviolet for dust mixtures appropriate to the diffuse interstellar medium in the 
Milky Way, the LMC, and the SMC. These scattering functions, which can be used for mod- 
eling reflection nebulae, are available at http://www.astro.princeton.edu/~draine . 
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2. The polarization of scattered light, as a function of scattering angle, is calculated at selected 
wavelengths. 

3. The calculated phase function 3>(0) for the WD01 model for Ry = 3.1 MW dust is compared 
to the widely-used Henyey-Greenstein phase function (po(0). At wavelengths 0.48/um < A < 
0.96/um, (f>o provides a good fit to <&, with r.m.s. relative error h re \ < 10%. In the ultraviolet, 
however, <pQ provides a poor fit, with r.m.s. absolute error 

/j abs > 50% at A < 2400 A. The 
Henyey-Greenstein phase function is not suitable for accurate modeling of reflection nebulae 
in the ultraviolet. 

4. A new phase function, 4> a <i(0) is proposed [see equation (5)]. This phase function provides 
a good approximation to our calculated phase functions for A > 0.27/um, with r.m.s. relative 
error h ie \ < 10%. At shorter wavelengths, however, $ becomes strongly forward-throwing, and 
neither the Henyey-Greeenstein phase function <f>$ nor the new phase function 4> a <i provide 
good approximations. 

5. There are significant discrepancies among observational determinations of the albedo and 
(cos 6) in the ultraviolet for the diffuse galactic light and dust in discrete clouds and reflection 
nebulae. One possible cause for these discrepancies may be reliance on the Henyey-Greenstein 
phase function in the radiative transfer models. It is recommended that the new phase 
function 4> a <i be used at wavelengths A > 0.27/um; at shorter wavelengths accurate radiative 
transfer models should use some other procedure to obtain <E>. 



I thank Karl Gordon, Aigen Li, Kalevi Mattila, Jonathan Tan, Adolf Witt, and Michael 
Woolf for valuable comments, the anonymous referee for helpful suggestions, and Robert Lupton 
for making available the SM software package. This work was supported in part by NSF grant 
AST-9988126. 



A. Phase Function Moments 



A phase function of the form 

M0) 



i 

47T 



i-g 2 



1 + a cos 2 9 



l + a(l + 2 5 2 )/3j {l+g 2 -2g cos Of/ 2 
has first and second moments given by 

l + a(3 + 2 5 2 )/5 



(cos 9) 
(cos 2 6) 



y l + a(l + 2 5 2 )/3 ' 
1 + 2g 2 + (3a/35) (7 + 20g 2 + 8g A ) 
3 + a(l + 2g 2 ) 



(Al) 

(A2) 
(A3) 
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B. Reproduce (cos 9) and (cos 2 9) 



If we specify (cos 9) and (cos 2 9), the adjustable parameters g and a are determined from 
equations (A2,A3). It can be shown that g satisfies the cubic equation 



g 3 - y (cos % 2 + 7(cos 2 0)9- I (cos 0) = . 



Let 



a ee I(cos 2 #>-|?(cos0> 2 



119 

6 ee ——(cos 6*) (cos I 
18 



If a 3 + 6 2 > 0, the solution is 



( a 3 + 6 2)V2_ 6 



1/3 



4913 / m3 7 / 
— — (cos6>) (cos I 

729 X 1 6 X 



( a 3 + 6 2 ) l/2 +6 j 1 /3 + l^ (co ^ ) 



whereas if o 3 + b 2 < 0, the appropriate root is 

g = 2|a| 1/2 cos(V>/3) + y cos6> , V = arccos (-&/l a | 3/2 ) 

The parameter a is then obtained from equation (A2): 

25((cos6>) - g) 



a 



3(3 + 2g 2 )g - 5(1 + 2 5 2 )(cos 9) 



(Bl) 

(B2) 
(B3) 

(B4) 

(B5) 

(B6) 



C. Fix a, Reproduce (cos 9} 
If we choose to fix the value of a and (cos 9), then g must satisfy 



g 3 -^(cos9)g 2 + 



3a + 5\ 5 , 

— S -^(3 + a)(c„ sS >=0 . 



For a / 0, equation (CI) has the real root 



9 



1/3 



(4 + ^) 1/2 -^] 1/3 + ^(cose) 



1 5 25 / „x2 

a Q - 2 + 6^-8T (c ° S ^ ' 

125 5 
- — (cos*> 3 + -(cos*> 



(CI) 

(C2) 
(C3) 
(C4) 



For a — ► (the Henyey-Greenstein phase function) we see immediately that equation (CI) requires 
9 = (cos 9). 
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